%% set parameter
steps = 200;
alpha = 0.0008; % learning rate
beta = 0.02;
stop = 0.0001;
rank = 600;
subRank=320;

%% run all matrix SGD

[ i ] = SGD( RateMatrix,rank,steps,alpha,beta,stop );

%% run the submatrix SGD
%the first submatrix
[Rec] = SGD( RecMatrix,rank,steps,alpha,beta,stop );
 [ RecoveryRate ] = RecoveryMatrix( Rec,indexLocation,n );